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How to use these programs 


Each program is arranged as follows: 


1. On the left of the page, explanatory information and the 
“execution sequence’, the sequence of keystrokes necessary for 
running the program. Results disptayed are printed in gold. 

2. Inthe first column on the right hand side of the page, the 
sequence of keystrokes which make up the program. 

3. In the second and third columns on the right hand side of the 
page, the program in check symbol and step number form 
(see section on checking the program). 


Notes 


1. Where a key has more than one function, the relevant function is 
printed as the keystroke in the first column 
cos 


e.g. the keystroke may appear as 8, cos or arccos. 


arccos v 
2. Thesymboi v within a program always refers to the key 
ChN/# 
3. The symbol # refers to 


4. The abbreviation gin is ‘go if neg’ and so refers to the key [1] 
go if neg 


Entering the program 


To enter a program into the calculator: 


1. Press [av] [ay] [0] [0] Display shows step programmed 


go to at 00 in check symbol form as 
described below. 


learn 


2. Press [av] No change in display. 


3. Press the sequence of keys for At each stage the step about to 
the program as shown in the __ be overwritten is displayed. 
first column of the program When the machine is first 


page. switched on every step is zero. 
4. Press Normal number display is 
resumed. 
5. Press [av] [av] [2] [o] [0] The step programmed at 00 will 
go to be displayed. 


Checking the program 
Each of the programs in the library is shown in check symbol form in 
the second column on the right-hand side of the page. 

Ss 


Press [av] repeatedly, and at each stage the check symbol will 
appear on the left of the display with the step number on the right. 
Ignore the four zeros in the display. 


e.g. A.0000 03 
check step 
symbol number 


After stepping through the program, press 


[av] [4¥] [0] [0] before execution. 


go to 


Finally, press and the program is ready for use. 


Correcting the program 


If the check symbol for a particular step number is not as indicated 
in the last two columns of the program page: 


Te Press [av] [ay] 


go to 
followed by the step number if the appropriate step number is 
not already displayed. 


learn 


2 Press [av] 


3. Enter the correct keystroke. The display will then show the next 
step in the program. If this is also incorrect, enter the correct 
keystroke. At each stage, the step about to be overwritten will be 
displayed. 

4. When correction has been completed, press( C/ce J. Any step 

‘ which has not been overwritten will not be atfected. 


5. Press [4¥] [av] [o] [o] 


go to 
Note 


To restore normal use of the calculator after entering or checking 
the program, press 


Running the program 


Press the sequence of keys as shown in the program library in the 
execution sequence. Results displayed are printed in gold. 


EXTENSION OF 
RANGE OF 
TRIGONOMETRIC 
FUNCTIONS 


to-—17<0<qT 


Sine of any angle: 


; ees _ 6 
sin we where t = tan 


Execution: 
6/RUN/sin@ 


For 6 in degrees, insert / v / D>R / at start of 
program. 





EXTENSION OF 
RANGE OF 
TRIGONOMETRIC 
FUNCTIONS 


to-17<0<T7 


Cosine of any angle 





cos 6 = 1+ where ‘Gini 
1 +t? 2 


Execution: 
@/ RUN /cos @ 





EXTENSION OF 
RANGE OF 
TRIGONOMETRIC 
FUNCTIONS 


to—17<0<7 


Tangent of any angle 


2t a 
150 where t = tan 





tan 6 = 


Execution: 
6 / RUN /tan 6 





EXTENSION OF 
RANGE OF 
TRIGONOMETRIC 
FUNCTIONS 


to—7<0<7 


3 . 6 
sin, cos and tan using t = tan> 


Execution: 
6/ RUN /sin @/RUN/cosé@ / RUN /tan 6 





SINE AND COSINE 
OF ANY ANGLE 


Sin: use program on right 


Execution: 
angle in degrees / RUN / sine 


For radians version of program, insert 
/ ¥/ RD / at beginning and omit / = / = / at end. 


Cos: either use program on right and execute by 
/ av / av /goto/0/4/ angle in degrees / 
RUN / cosine 

or omit first four keystrokes of program 


on right and fill the empty spaces at the 
end with repeated / = / and execute by 


angle in degrees / RUN / cosine 


For radians version of program, insert / ¥ / R>D / 
at the beginning. 


Note: E can appear if reduced angles > 1-57 
radians. 





TANGENT OF ANY 
ANGLE 


Execution: 
angle in degrees / RUN / tangent 


Note: E can appear if reduced angle > 1-57 
radians. 





HYPERBOLIC 
FUNCTIONS 


If all the hyperbolic functions are likely to be 
required, use the ‘gudermannian’ program on 
page 21 . For the individual functions, the 
following can be used: 


Sinh x 


Execution: 
x / RUN / sinh x 


Range: 
—227:95 <x < 230-25 
(out-of-range may give wrong result without E) 
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HYPERBOLIC 
FUNCTIONS 


Cosh x 


Execution: 
x / RUN /cosh x 


Range: 
—227:95 <x < 230-25 
(out-of-range may give wrong result without E) 





HYPERBOLIC 
FUNCTIONS 


Tanh x 


Execution: 
x / RUN / tanh x 


Range: 


|x| < 113-97 
(out-of-range may give wrong result without E) 





HYPERBOLIC 
FUNCTIONS 


Cosech x 


Execution: 
x / RUN / cosech x 


Range: 
1:0017 x 10% < |x| < 227-95 
(|x| > 227-95 may give wrong result without E) 


>2] | |m[o[o]o/o/ 7] =>) 
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HYPERBOLIC 
FUNCTIONS 


Coth x 


Execution: 
x / RUN /coth x 


Range: 


1:0016 x 10% < |x| < 113-97 
(out-of-range may give wrong result without E) 
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HYPERBOLIC 
FUNCTIONS 


All the hyperbolic functions 


Execution: 


x / RUN /sinh x / RUN / cosech x / RUN / 
cosh x / RUN /sech x / RUN / tanh x / RUN / 
coth x / 


Range: 
10017 x 10% < |x| < 7-8566 
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HYPERBOLIC 
FUNCTIONS 


The gudermannian program 


Enables all the hyperbolic functions to be 
calculated with suitable execution sequences. 
Formulae: 

gdx = 2 arctan tanh 
sinh x = tan gdx 
cosech x = cot gdx 
cosh x = sec gdx 

sech x = cos gdx 

tanh x = sin gdx 


coth-x = cosec gdx 


Execution: 

x / RUN / gdx/ av /sin/ tanh x /+/=/cothx 
x / RUN / gdx/ av /cos/sech x /+/=/cosh x 
x / RUN / gdx/ aw /tan/sinh x /+/=/cosech x 
This program can be used inside parentheses and 
does not affect memory. 

Accuracy is less than that of individual 
hyperbolic function programs. 

Range: 

|x| < 227-95 for gdx 

|x| < 7-8566 for hyperbolic functions 
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INVERSE HYPERBOLIC z 
FUNCTIONS F 
3 04) 
_1 05) 
All the inverse hyperbolic functions can be 108 
obtained from the following program. 





Execution: 
av / av /goto/1/2/x/ RUN* /sinh |x 


or 
x / RUN /+cosh'x / RUN / —cosh' x 

or 

av / av /goto/2/0/x/RUN/tanh 'x 

or 

av / av /goto/1/3/x /RUN* /cosech'x 


or 

av / av /goto/0/1/x/RUN/+sech'x / 
RUN /—sech * x 

or 


av / av /goto/1/9/x/ RUN/coth 'x 


* For negative x press / RUN / a second time 
when evaluating sinh’ x and cosech™' x to get 
the correct answer. 


Range: 

sinh’ x 10°49 < |x| < 577-35 

cosh 'x 1<x <3162-2 No E if x —ve 
tanh x —0-99999 < x <0-99999 
cosech?x 0-001732 < |x| < 10*° 


sech"' x 3-162278 x 104 <x <1 
No E if x —ve 


coth'x 1:0001 < |x| < 1099 


"stop | 0 [32 
P [F338 
=| fea 
stop | 0 [355 
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MODULO ARITHMETIC 
(‘Clock Arithmetic’) 


Base 7 is used as an example. 


The program completes a calculation and works 
out the remainder when the result is divided 

by 7. Neither the brackets nor the memory are 
used, so that the operation of / RUN /is 
exactly that of /=/. 


For other bases, insert the base at steps 03, 11 
and 14. Change the address at steps 18 and 19 
to 14 if a two digit base is used, 16 for a three 
digit base, etc. 


Execution may take a long time if very large 
numbers are involved. 


Example: 
/3/X/5/RUN/1/+/8/ RUN / 2 / ete. 
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PRIME 
FACTORISATION 


To find the prime factors of a number N. 


Pre-execution: 
2/ av /sto/ av / av /goto/0/0/ “ce / 


Execution: 


N/ RUN /a,/RUN/N,/ RUN/2,/ RUN / 
No/+++/a,/ RUN /1 


where 

a1, 42, a3, °° * a, are the prime factors of N 
and 

N,, N2,-+- are the residues defined by 


ne eh gg a eN No ete. 


a, fire ax! scabs 6 





Pressing / RUN / after 1 has been displayed will 
cause the machine to go into an infinite loop. 


Warning: Long execution times are possible for 


large values of N or for numbers with large 
prime factors. 
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PRIME 
NUMBER TESTING 


To find whether a number n is prime, choose 
any integer m >/n. 

Then use the execution sequence: 

n/ RUN /m/ RUN / 

The result will be the largest number less than 
or equal to m which divides n. If the result is 1 
then n is prime. 

To test another number, pre-execute with: 

/ ay / av /goto/0/0/ 


Note: Long execution times are possible for 
large numbers. 








FACTORIALS 


Execution: 
n/RUN/n! 


Restriction: 
1<n<69 


Note: The program may be used within pws 3H 
brackets. It does, however, use the memory. 18 
Thus, to calculate A 

151 hand 

6! 10! = | 


@: 


a possible execution sequence is: 


15/RUN/+/ av /(/10/RUN/ av /) /+/ 
av /(/6/RUN/ av /)/=/ answer 
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FACTORIALS 
OF LARGE NUMBERS 


This program calculates In (n!) for n greater 
than about 25. 


Reasonably accurate results are given for 
n greater than 10. 


(The program uses, Stirling’s approximation, 


n! =,/2mne"n") 
Execution: 
n/ RUN / In (n!) 


_stop_| 0 [20] 
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THE GAMMA 
AND PI FUNCTIONS 


IT (n+ 1) =27(n) =n! when n is an integer 


I (x + 1) = x I(x) for x >0 

I’ (0) is undefined m (0) =I (1) =1 
lr (4%) =r T (%) = h/t 
I (1)=1 m (1) =I (2) =1 


By interpolation: 

lr (n+6)#(n+%d—%)Pr(n) OK< 

‘ Pe 5 1 (r) 

“mt (6) = (n+ %6 — %) H 5 OS 

ate) na (nthd—%) nor 
6 r=1(r + 6) 


1 


6< 
6<1 


I (5) = 





O0<6<1 


n should be suitably large for the accuracy 
required. 


n= 20 gives high accuracy 

n=5 gives reasonable accuracy for most 
purposes 

e.g. T (2) = %./m = 0-8862269 

n=5_ gives m (%) =0-885547 

n= 20 gives m (%) =0-8861174 


Execution: 


av / av /goto/0/0/ 


n/RUN/&/RUN/n—1/RUN/n—2/ 
RUN /---/2/RUN/1/RUN/ «(5) / av / 
av /goto/3/2/RUN/! (5) 
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FIBONACCI NUMBERS 


Each number in the sequence is the sum of the 
previous two. 


Execution: 
Cce / RUN /F,/RUN/F,/ RUN /--- 





NUMBER 
BASE CONVERSIONS 


Decimal to binary (fractions) 


Given a decimal x, O< xX <1, this program 
calculates the binary expansion of x to any 
number of places. 

Suppose x = 0-d,d, . (binary) 


Execution: 
x / RUN/d,/RUN/d,/RUN/d,/--- 


To calculate the expansion of another decimal 
y, press 


[ce /Cce/ av / av /goto/0/0/y/RUN/--- 


ete; 


Notes: 


1. To convert decimal integers to binary use 
the program on page 31. 


2. No program for converting decimal 


fractions to bases other than 2 is provided. 
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NUMBER 
BASE CONVERSIONS 


Decimal integer to base m 


This program expresses any integer in any base. 


Suppose x = a, -:: a, in base m. 


Execution: 


m/RUN/x/RUN/a,/RUN/a,_ ,/---/a,/ 
RUN /m 


Note that the digits are produced in reverse 
order and that the machine tells you that all the 
digits have been shown by displaying the base m. 


The sequence can be repeated for a new x and/or 
m. If the same m is required there is no need to 
re-enter it because it is already in the display. 


Note: To convert decimal fractions to base 2, 
use the program on page 30. 
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NUMBER 
BASE CONVERSIONS 


Binary to decimal (integers) 


Binary isa,--- a, 


Execution: 
a,/ RUN /a,_,/RUN/--+/a,/=/ answer 
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NUMBER 
BASE CONVERSIONS 


Binary fraction to decimal 


If number is: 
0: bb, -:* by 


Execution: 
RUN /b,/ RUN /b,/--+/b, / RUN / eoswer 


At each stage the answer so far is displayed. 


Fraction base m to decimal 


Exactly the same except / 2 / at step 10 is 
replaced by the appropriate base. 





NUMBER 
BASE CONVERSIONS 


Binary to decimal (integers, fractions or 
mixed numbers) 


Binary isa, 3..a> * bys«- Dp, 


Execution: 


Gce / RUN /a,/ RUN/a,_,/RUN/---/ 
RUN /a,/—/ RUN /b,/ RUN /b,/---/b,, / 
RUN / answer 


Notes: 


1. The/—/ corresponding to the ‘decimal’ 
point must be entered even if the number 
is an integer. 


2. The correct answer will be given if: 
a, = 1 
n21 
a, =1or0 


orjust-b,--- (- entered as —) 


To re-use: 
Cce /Cce/ av / av /goto/0/0 


34 









MEx 5 (26 
|) {6 [27] 
| vA [28] 
| 1 {1 [30] 
| 6 {6 [31] 
| || 
| | (33) 
| | 84 
Lt Bs 










NUMBER 
BASE CONVERSIONS 


Base m to decimal (integers) 


Number isa, a,_4,°** a 


Execution: 
m/RUN/a,/RUN/a,_,/RUN/---/a,/ 
= / answer 

To re-use with same m: 

Cce / RUN /a’ --- 
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NUMBER 
BASE CONVERSIONS 


Base m to decimal (integers, fractions or 
mixed numbers) 


Number isa... 5a, *'by... Dp, 


Example: 


m=7 


Execution: 


Cce / RUN /a,/ RUN/a,_,/---/RUN/a,/ 
—/RUN/b,/ RUN/b,/---/RUN/b,/ 
RUN / answer 


Notes: 


1. Insert value of m at 02 and 25. 

2. If two digit base is used, insert at 02, 03, 
move the next 22 steps down one, insert 
the base again at 26, 27, and substitute 
/1/9/for/1/8/ in the last two steps. 
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SERIES 





Natural numbers 





PF PEC EER ARR RRRRREERS 


n(n + 1) 


4 
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-'+n)= 
n/ RUN /sum 


(1+2+ 
Execution: 
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SERIES 


Squares of natural numbers 


(144494 0-412) = n(n 1)(2n+ 1) 


Execution: 
n/RUN/sum 
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SERIES 


Cubes of natural numbers 


(148 +27+++ +n) =2n2(n +1)? 


Execution: 
n/ RUN / sum 
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ARITHMETIC SERIES 


First term =a 
Common difference = d 
N terms 


(N-— ue) 


= + 
sum N( <2 5 


Execution: 
a/RUN/N/RUN/d/RUN/ so 


40 










+ 





ale 
oD 
TNO OM 


| 
BRRBSLS 






ne 
ete 
g 


A 
Ce 






4 
x ° 

TC 
mk | ek 
SSso8 








n 

or 

o 
<Bea 
Ea al 
— 

5 


a 
° 
os 

3 
Ha 







NNN 
yO - oO 


NM NM NY NY 
BaRS 





NO 
~N 






NO 
© 






WwW WwW 
WON — 


= 
| | 85) 


ARITHMETIC SERIES 


‘First term =a 
Last term = | 
N terms 

a N(a + |) 
Z 


Execution: 
a/RUN/I/RUN/N/ RUN /sum 





GEOMETRIC SERIES 


a(1 —r") 


=at pene t+ N-1— 
S=atar ar (=n 


First term =a 
Common ratio = r 


N terms 


Restrictions: 
r>0, r#¥1 


Execution: 
a/RUN/r/RUN/N/RUN/sum 
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INFINITE 
GEOMETRIC SERIES 


S=atartar?+--- = 





1-—r 
Restriction: 
Ir}<1 


Execution: 
a/RUN/r/RUN/sum 
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SUMMING SERIES 
IN GENERAL 


N 
z a(n), some function a. 


Examples: 
1; +4449 4--- +N? a(n) =n? 
& 1) 4. ss gf Wee 
(1+7)+ (8+3)+ +(N Re) 
a(n) = n? pe etc. 
n 


Write a program segment which evaluates a(n) 

when n is in memory; parentheses may not be 

used. The segment may be up to 15 steps long, 

any final / = / stop / being omitted. Fill up any 

unused steps with /=/---/=/. 

Examples for above: 

ter on? rel / X / 

a net write as (n° + 1) +n? 
rel /X/X/X/rel/+/#/1//rel/ rel / 
=/=/=/ 


Then use the program as shown. 


Pre-execution: 


Clear memory with Cre / av / sto / 


Execution: 
N/ RUN / a(1) + a(2) +--+ - + a(n) 





SERIES 


a(x,) +a(x2) +--+ +a(x,) 

Write a program segment to evaluate a(x;) 
without using parentheses; the memory may 
be used. 


Then use the following program: 
/(/stop/---seg---/)/+/¥/goto/0/0/ 


Execution: 
RUN /x,/ RUN /x2/---/x,/ RUN / sum 


At each step the sum so far is displayed. 
Example: 
To find ~ tan (x2 +1) 


1 x3 +1 
Express x? +— as ——— 
x x 


Program segment is then: 
/sto/X/X/rel/+/#/1/+/rel/=/tan/ 
and so program is as shown. 

The segment may be up to 32 steps long, by 


omitting / ¥/ goto /0/0/ at the end and 
filling any empty steps with / = /. 
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HARMONIC ADDITION 


Resistors in parallel, capacitors in series, lenses 
in series, etc. 


1 

—= i + . + — 
X = 2X, x. 
Execution: 


x, / RUN / x2/ RUN / +++ /x,/ RUN / x 


At each step the harmonic sum so far is 
displayed. 





PYTHAGOREAN 
ADDITION 


Geometry, electricity 


x= Vx? Fee Xe 
Execution: 
X,/ RUN /x./ RUN /---/x,/ RUN / « 


At each step the intermediate result 
xj +-+--+x? is displayed. 
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ARITHMETIC MEAN 


Pre-execution: 
Gece /Cce/ av / av /goto/0/0 


Execution: 


X,/ RUN /x/ RUN /---/x,/ RUN / 
arithmetic mean 


At each stage the arithmetic mean so far is 
displayed. 
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GEOMETRIC MEAN 


Pre-execution: 
Cce/oce/ av / av /goto/0/0 


Execution: 


x,/RUN/x,/ RUN /---/x,/ RUN / 
geometric mean 


At each stage the geometric mean so far is 
displayed. 
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HARMONIC MEAN 


BEN diel fh) 
H n\x, Xn 
Pre-execution: 

Gce /Cce/ av / av /goto/0/0 


Execution: 


x,/ RUN /x,/RUN/---/x,/ RUN / 
harmonic mean 


At each stage the harmonic mean so far is 
displayed. 
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ROOT MEAN SQUARE 


JV (x2 +++ + + x2) 


n 


R = 
Pre-execution: 
Cce/Sce/ av / av /goto/0/0 


Execution: 


x,/ RUN /x2/---/x,/ RUN / 
root-mean-square 


At each stage the r.m.s. so far is displayed. 
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QUADRATIC 
EQUATIONS 


ax? + bx +c=0 
Roots x,, X2 if real 
R + il if complex 


Execution: 

x, /RUN/»,./ RUN / 
RUN / Cre / Cre / if roots 
are real 

1° /See /RUN/ 8 / 

if roots are complex 


a/RUN/b/ 
RUN /c/ RUN / 


* error symbol displayed 


After the sequence a/ RUN /b/ RUN /c/ 
RUN / the display shows e/ther (if the roots are 
real) the larger real root with no error indication 
or (if the roots are complex) the imaginary part 
and the error symbol. Continue with the 
appropriate execution sequence. 


The error symbol will tell you whether the roots 
are complex. The sequence / RUN / RUN / Cre / 
shown above after (x) is necessary before 
entering a new equation to be solved. 









‘stop 026 
| ret_| 5 [271 
_rel_| 5 |29| 
| =] = [304 
stop_0 31. 
vx | 1/32) 
He {Jas 
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stop 0 36, 
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CUBIC EQUATIONS 
by an iterative method 


ax? + bx2+cx+d=0 


Formula: 


Ot. 8 2axz +bxz—d (based on Newton- 
k +1" 3ax2 + 2bx, +c Raphson method) 


(Fill in your own values of 2a, b, d, etc.; 
if any of these are negative change the + or — 
preceding them to — or +) 


Execution: 


. d 
Choose any starting value x,, say Ee 


x, /RUN/x,/RUN/x./--- 


If the sequence converges, the limit will solve 
the equation. 


If the sequence does not converge, try a new 
starting value. 


The sequence will usually converge to the root 
closest to the starting value and so by trying 
different starting values all the roots should be 
obtained. 


* where a, a, is the two digit number 3a; 
if 3a << 10 then enter a, = 0 and a, as the value 
of 3a. Similarly b,b, is 2b. 
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POLYNOMIALS 


To evaluate 


a,x"ta,_4x"7! +---+a, = p(x) 


Execution: 


x /RUN/a,/RUN/a,_,/:--/a,/RUN/a,/ 
= / result 


To use again: (with different x) 
ay / av /goto/0/0/ before execution 


Notes: 


1. The individual results after each / RUN / 
are the coefficients of the polynomial q(x) 
where q(t) = p(t) / (t — x). 

2. If p(x) =0, x is a root and q(x) is the 
quotient polynomial which can be solved 
for other roots of p(x). 








POLYNOMIALS 


To write a program to evaluate the same 
polynomial repeatedly 


Example: 
p(x) = 5x4 + 8x3 — 3x2+ 4-2x +1 


Method: 
Express as | [(5x +8)x —3]x+4-2]x +1 


Execution: 
x / RUN / p(x) /y/RUN/p(y) ++: ete. 


Note: \f a coefficient is zero omit it together 
with the — or + sign preceding it. If the leading 
coefficient is 1, it may be omitted together 
with the multiplication sign which precedes it. 
See over for example. 
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POLYNOMIALS 


Example: 
To calculate x? + 2x?+3 


first coefficient = 1, 
so omitted. 


coefficient of x =0, 
so omitted. 











DIVISION OF 
A POLYNOMIAL BY 
A QUADRATIC 


Division of the polynomial 
p(x) =a,x"+ a, yx" +--+ +axta, 
by the quadratic divisor 
d(x) = x2+mx+n 
gives the quotient polynomial 
q(x) = b, 2x"? +b, gx"9+---+b x +b, 
with remainder 
r(x) =c,x +c, 
Pre-execution: 

av / av /goto/0/0/ceE/ av /sto / 
Execution: 
RUN /n/RUN/m/ RUN /a,/ RUN / 5. . 
RUN /n/ RUN /m/ RUN /a,_,/RUN/b._, 
RUN /n/ RUN /m/RUN/a,/RUN/b, 
RUN /n/RUN/m/ RUN /a,/RUN/c, 
RUN /n/RUN/m/RUN/a,/ RUN / RUN / 
m/RUN/1/RUN/c, 
/ RUN / RUN / completes execution 
Results may be tabulated as below: e.g. to divide 


x® — 4x5 + 31x4 — 96x? + 415x? — 652x + 1105 
by-x? + 2x + 3: 


r n m a, b,_o 

6 3 2 1 1 

5 —4 —6 

4 31 40 

3 —96 —158 

2 415 611 

1 —652 —1400=c, 
0 1105 —728 =c, 
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SOLVING 
A POLYNOMIAL 


This is an iterative method to find a quadratic 
factor of a polynomial. When the polynomial 
has been reduced to quadratic factors, these can 
be solved to give the real or complex roots of 
the original polynomial. 
Stage 1: 
Choose a starting quadratic divisor 
d(x) = x? +mx+n_ (say) 
Divide p(x) by d(x) to give a quotient q(x) and 
remainder r(x) = rx +s 
Stage 2: 
Divide q(x) again by d(x) to give a new quotient 
q'(x) and remainder r’(x) = tx +u 
Stage 3: 
Find the coefficients m’ and n’ of the next 
iterate of the quadratic divisor using this program 
Execution: 
u/ RUN /t/ RUN/m/RUN/n/RUN/1D 
p/RUN/t/RUN/u/RUN/s/RUN/r/ 
RUN /t/RUN/—/+/n/=/9 

av / av /goto/2/5/r/RUN/u/RUN/ 
n/X/s/RUN/t/RUN/+/m/=/m 


D = u2 + nt? — mut 


d= @ + Ent 
D 
; rt +s(mt — u) 
n =n - —_ 
D 


Re-enter the quadratic divisor program and 
iterate again with the new values of m’ and n’. 
Repeat stages 1—3 until the values of m and n 
converge. 
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NUMERICAL 
INTEGRATION 


Triangular interpolation 


' 


1 


2 
Execution: 


oe 2Vn-41 + Yn) 


hly, tf 2y; = 2y2 io 





n/ RUN /y,/ RUN /y,/ RUN / 
/RUN/y,/RUN/h/ RUN/|! 


° 
pul io 
eee 7 


£ 
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NUMERICAL 
INTEGRATION 


Simpson’‘s Rule 


1 
I= Zhlyo + Ay, + 2y2 + 4y3 + FAV get Val 
(n must be even) 


Execution: 


n/—/1/=/RUN/y,/RUN/---/y,/RUN/ 
h/ RUN / | 





| 0 | 0/25) 
stop 0 27 
Lh] 6128) 






== faa 
“stop 0 35: 
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INTEGRATION 


Weddle Formula 


NUMERICAL 


(Yo + 5y, + yo + Gy3 + Ya + 5ys5 + Ye) 


_ 3h 
10 


Integral 


Execution: 


/ RUN/y2/ RUN/y3/ RUN /y,/ 


/ RUN/yg/ RUN /h, / 


Yo/ RUN /y, 
RUN /ys 
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COMPLEX NUMBERS 


z=xtiy 


To find magnitude and argument. 


Execution: 


If y = 0, then z = |x| and arg z = (O if x >0, 
t if x <0) 


Otherwise, x / RUN /y/ RUN/ |z\/ RUN / arg z 
To find x and y given arg z and |z| 

(—1 Sarg z <7) 

If arg z is O, then x = |z| and y=0 


lf arg zis 7, then x = —|z| and y=0 


Otherwise use polar-cartesian program, 
execution as follows: 


\z| / RUN /argz/ RUN / x/ RUN/ y 





DETERMINANTS 





a, b, 

=a,b,.—a5b 
a> | 412 2s 
Execution: 


a, / RUN /b,/ RUN /a,/ RUN /b,/ RUN / 
det 
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MATRIX 
MANIPULATION 


1. Matrix multiplication (steps 0O0O—11) 
AB=C . 


C; = 


Execution: 

aj, / RUN /b,; / RUN / 

aig / RUN /b,/ RUN /--- 

ain / RUN /b,, / RUN / ©, 

To restore zero total for next calculation, 
press CE, 

Error correction: 

Foray: CeE/+/a, 

Forb,j: av /)/SKeeE/+/ av /(/ by; 


2. Back substitution (steps 00—21) 
(for AX = B where A is upper triangular) 


ij = 
Pre-execution: 
av / av /goto/0/0/ for each x; 
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MATRIX 
MANIPULATION 


Execution: 
Xnij / RUN / ain vf RUN /- = /Xiaay if RUN / aiisd i RUN / Dain Xk; 
av / av /goto/1/2/ RUN /b,/ RUN / a; / RUN / x. 


_ Error correction: 

For xyj: SCE /+/ xX; 

Fora,: / av /)/S%ce/+/ av /(/a, 
Forby: Se /—/b; 

Foraj: Sce/+/a, 


3. Adding a multiple of row i to row j in the augmented matrix (A/B) 
(steps 16—30) 


Aik = Aj + Max, dix = Dix + Mii dix 
ae 
ji 


where m,, = —-— 
aii 


Pre-execution (each m;;): 
av/ av /goto/1/6/CKeE 


Execution: 

aj / RUN / a, / RUN / ™,, error correction: re-run from 1¢€ 
ax / RUN / a, / RUN / a), for each k 

b;. / RUN / bi. / RUN / bi, for each k 


Note: \f mj, is known pre-execution canbe av / av /goto/1/9/ ct 
and first part of execution m,; / RUN / 7, 
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EQUATION SOLVING 


The secant method 


In this variant of the Newton-Raphson method 
for solving the equation f(x) = 0, instead of 
computing the derivative f’(x) at each stage, an 
approximation to f’(x) at a point in the vicinity 
of a root x, is used. 


Stage 7: 


Write a program segment to compute f(x) when 
x is in memory, taking up no more than 27 steps 
excluding the final / stop /. Enter the program 
starting at step 01, ending with the sequence 

/ stop / ¥/ goto /0/0/. 


x / RUN /¢ (x) 


Evaluate f(x) for a range of values in which a 
root is likely to occur. If f(x,) and f(x.) have 
opposite signs, there is a root between x, and x. 


Execution: 


Stage 2: 


Calculate an approximation to the derivative of 
f(x) as follows: 

f(x2) /—/ (x4) /+/ av /(/x,/—/x5/ av / 
)/=/k =—'(x,) 


Stage 3: 


The iteration formula for the secant method is 


x + 


where K is a constant approximately equal to the 
derivative of f(x) at the root. K may be chosen 
to be equal to k, or may be an integer or a 
number with fewer digits than k, in which case it 
should be numerically larger than k. 


x’ = 


Note: \f the program segment in Stage 1 took 27 
steps, there is room for only one digit for K in 
the following program. (contd. over) 
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Starting at the final / stop / step, press / av / LEARN / and enter 
the sequence: 


/+/#/K/+/rcl/=/stop/ for K positive, or 
/-—/*+/#/K/+/rel/=/stop/ for negative K 


The sequence / vw / goto /0 /0/ or /=/ steps may be added at the 
end. : 


Execution: x/RUN/ 
/RUN/* -- 


Repeat until successive values are equal. If convergence is slow, 
‘decrease K. If the results diverge, increase K. 


If k is a small fraction, the / + / step may be replaced by a/ X / step 
. and K taken as the reciprocal of k. 


See below for example. 
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Example: 
To solve cos x = x 
f(x) = cos x — x 





us 
Take x, =9" Xo = 0. 
T 
1+>= 
Then {%o) — fly) 1 _ 2 2, 
X1— Xo 1 
2 


Program segment is / cos / —/ rel 
Guess 1 as starting solution 
Execution: 


1/ RUN /0:770223 
/ RUN /0-7440342 
/ RUN / 0-7399375 
/ RUN /0-7392705 
/ RUN / 0-7391738 
/ RUN /0:7391592 
/ RUN /0-7391519 
/ RUN / 0-7391483 
/ RUN /0-7391465 
/ RUN /0-7391456 
/ RUN /0-7391451 
/ RUN /0-7391449 
/-RUN /0-7391448 
/ RUN /0-7391447 
/ RUN /0-7391447 


So result is 0-7391447 








CIRCLES 


Circumference and area 


N ~ N 
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Execution: 
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CIRCLES 


Radius of circle from area 


a|> 


Execution: 
A/RUN/+r 








N 


NN NIN NM NM NM 
BL SRRBS8S 











CIRCLES 


Area of segment: 


ea? 


Area of segment if h and r are given: 


r2 
Area => (a — sin a) 


a r—h 
where cos = = ——— 
2 r 


Note: the angle a is calculated internally and is 
not required to be input. 


Execution: 
r/RUN/h/ RUN / arco 
Note: limited range, a < 1-57 radians 
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CIRCLES 


Length of chord 


a? 


|= 2 /2hr — h? 


Execution: 
h/RUN/r/ RUN /\enoth 











CIRCLES 


Area of circular annulus 


tee 


Area = 1(R2 — r?) 


Execution: 
R/RUN/r/ RUN / area 
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TRIANGLES 


To find area, given base and height 


Execution: 
b/RUN/h/ RUN / area 











TRIANGLES 


To find area, given all three sides 


Area =_ /s(s — a)(s — b)(s — c) 


apres) 


where s = ( 2 


Execution: 


a/RUN/b/RUN/c/RUN/b/ RUN /a/ 
RUN / area 
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TRIANGLES 


Finding a side, given two angles and a side 





_bsina 





sin B 


Execution: 
a° / RUN / 6° / RUN/b/ RUN /a 








TRIANGLES 


Length of third side from two sides and 


included angle 
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b? + c? — 2bc cosa 


a= 


Execution: 
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b/RUN/c/RUN/a@°/ RUN / 2 
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TRIANGLES 


Finding an angle, given three sides 


| ax 


b? + c? —a? 


cos a= Sinn 


Execution: 
a/RUN/b/RUN/c/RUN/« 
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Length of medians, given lengths of sides 


TRIANGLES 





See seer eee RANA AVAL SIS 8 





2(b? + c?) — a? 





2 
Execution: 
b/RUN/c/RUN/a/RUN/S 
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80 
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RIGHT ANGLED 
TRIANGLES 
Length of hypotenuse from other two sides 
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Execution: 








a/RUN/b/RUN/) 
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TRIANGLES 
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h2 — a2 
a/RUN/h/ RUN/5 


Execution: 
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PARALLELOGRAMS 


Area = ab sin 0 


Execution: 
a/RUN/b/RUN/6° / RUN / area 


For @ in radians, insert / ¥ / R>D / between 
steps 04 and 05. 


goto |2 |17 
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Anr? 
mr? 


+1 


radius / RUN / surface area / RUN / volume 


SPHERES 
Surface area and volume 
Execution: 
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SPHERES 


Radius from volume 


Execution: 
V/RUN/+r 











CYLINDERS 
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= rh 
Area of curved surface 


Volume 


. 


2mrh 


Total surface area = 2zr(r + h) 





Execution: 


r/RUN/h/ RUN / volume/ RUN / area of 
curved surface / RUN / total surface area 
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HOLLOW 
CYLINDRICAL TUBE 









Area of curved surface = 27h(R + r) 


Volume = zh(R?2 — r?) 





| 6/6421) 
Execution: oe A 
R/RUN/r/RUN/h/ RUN / area of curved stop O 23 
surface / RUN / volume rel = 5 24 

om. EP 26) 
stop 0 26 
| v_ |A(27| 
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RIGHT CIRCULAR 
CONE | 


2 

Volume = ae 

Curved surface area = mr. /r2 + h? 
Total surface area = mr (r +V/1? +h?) 


Execution: 


h / RUN /r/ RUN /area of curved surface / 
av /rcl /area of base / RUN /total surface 
area / RUN / volume 


88 







-* 
















wn 
+t +x 3 8 
oS 
Ssssesss 
~S 


ar 
ene 













> 
a 


aes 
> [eo] eof «| [ey 
Bass 


— 
N 







a> owas 
[3/3 /3/3/8/8 [Slee 
N > 





n 
+ 
9° 
To 


Doo 
288 








p 


n 
ot 


RECTANGULAR 
PARALLELEPIPED 





Diagonal: 


d=, /a?+b?+c? 


Execution: 
a/RUN/b/RUN/c/RUN/ © 





RECTANGULAR 
PARALLELEPIPED 


Surface area 


A = 2(ab + ac + be) 


Execution: 
a /RUN/b/RUN/c/ RUN / area 
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DISTANCE BETWEEN 
TWO POINTS IN SPACE 


d= /(x,— Xz)? + (y, — yo)? + (z, — 22)? 
points are (x,,y;,Z,) and (x2,y2,Z2) 
Execution: 


x,/ RUN /x,/ RUN /y,/ RUN /y2/ RUN / 
z,/RUN/z,/RUN/¢ 
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COORDINATE C= Te Too 
CONVERSION £8 ot 


Cartesian to polar 


Restriction: y #0 

If y=0,r=I|x|_ 
and 6 =Oifx20 

| nifx <0 


Execution: 
x /RUN/y/RUN/: /RUN/~ 
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HAVERSINE AND 
INVERSE HAVERSINE, 
VERSINE 

AND SUVERSINE 


Haversine: 


pre-execution: av / av /goto/0/0/ 


Execution: 
6° / RUN / hav 0 
/+/=/vers0/—/+/2/=/ suvers 0 


Inverse haversine: 


pre-execution: av / av /goto/1/4/ 


Execution: 

hav 6 / RUN / 0° 
vers@/+/2/=/RUN/0- 

suvers 0 /—/+/2/+/2/=/RUN/0- 
Range 0< 6° < 180 


For vers 6 see post and pre-execution. 





AREA OFA 
SPHERICAL 
TRIANGLE 





Area = (A+B +C — 2z)R2 
A, B,C in degrees 


Execution: 
A/RUN/B/RUN/C/RUN/R/RUN/~.- 
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SPHERICAL 
TRIANGLES: SINE RULE 


sinA_sinB_ sinC 


sina  sinb - sinc 


Angles in degrees 





Execution: 


a/RUN/A/RUN/B/RUN/) /C/ 
RUN /c 

or 
A/RUN/a/RUN/b/RUN/® /c/ 
RUN /C 


Note: \|f a result of O appears, the final arcsin 
had an out-of-range argument and the result is 
impossible for the particular angles given, or 
else very close to 90°. 

For angle A > 90°, compute using 

180 /—/A/=/ete. 


Special execution: navigation 


To find course from place 2 to place 1 


sinc < sin (E, : E.) cos No 
sind 


Execution: 


E,/—/€,/RUN/d/RUN/90/—/N,/=/ 
RUN /C 
where E, = easterly longitude of place 1 

E, = easterly longitude of place 2 

Nz = north latitude of place 2 


d = angular distance between places 
1 and 2 


(For westerly longitudes or south latitudes, 
change sign of angle.) 





SPHERICAL 
TRIANGLES: 


Cosine Rule 


om 
(D 
ey 


cos a= cos b (cos c+ sinc tan b cos A) 


Execution: 
c/RUN/A/RUN/b/RUN/b/RUN/ 2 
Navigation 

To find great circle distance between places 
1 and 2 


1. Latitude N, longitude E, (—ve if W) 
2. Latitude N, longitude E, (—ve if W) 


Execution: 


90 /—/N,/=/RUN/E,/—/E,/RUN/90/ 

—/N,/RUN/90/—/N,/ RUN / a (degrees) 
X /111-19 /=/ distance in km 

or X/69-41 /=/ distance in miles 


For angles greater than 90° use appropriate 
reductions to first quadrant. 
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= |= bo] 
yA Bt 
arccos 8 32 
Py [A faa) 
Fed 6 34) 
‘stop /0 35) 








SPHERICAL 
TRIANGLES 


The Cosine Rule — to find an angle or side given 
three sides or angles 





cos a —cosbcosc 


cosA= : : 
sinb sinc 


cos A+ cosB cos C 


cos a= = : 
sin B sin C 


Execution: 
Angles in radians 
c/RUN/b/RUN/—/RUN/a/RUN/A 
C/RUN/B/RUN/RUN/A/RUN/a 
For angles in degrees use 

av / av /D>R/ after each angle. 






__sto | 2/00 
01 


|_sin | 7 
per 1S) 
_MEx 5 03 
cos 8 
ee 
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| stop | 0/07, 
sin 7 08 
ea) 
LP 2¥_ PAO 
_MEx 5 11 
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SPHERICAL 
TRIANGLES: 


Half-angle tangent formula 


A sin (s — b) sin (s —c) 
aS /- 
2 sin (s — a) sins 

ab +c 
where s =~ 


tan2=- /eos(*—S) cos (S — A) 
2 cos (S — B) cos (S — C) 
wheres = Se 
2 
Execution: 
Angles in radians 
a/+/b/+/c/RUN/b/RUN/a/RUN/ 
A 
5 /=/4 


For the cosine version change all / sin / steps 
to / cos /. 
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“stop 035 
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SPHERICAL Ee 

TRIANGLES T= [6 joa] 

Solving the PZX triangle - on ate 

| xX | + [05) 

| X_ | [06] 

| ( [6 [07] 

stop 0 08 

sto 2 09 

_sin 7 10 

|) [6 [11] 

Ke 

| | ( [6 [13] 

hav ZX = hav (PX ~ PZ) + sin PX sin PZ hav ZZPX ea 

= hav (L ~ D) + cos L cos D hav ZZPX | ¥ [A /16| 

“(for the second formula use cos at steps 10 and /MEx | 5 [17, 

27 instead of sin) | + [G ]18) 

ZX is the calculated zenith distahce (CZD) se 

Enter south latitudes as —ve tae he at 

Execution: 
Angles in radians — 2 

| = [= [24] 

LZPX / RUN / PX / RUN /PZ/ RUN /+/=/ 2% | wv [A [25] 

Angles in degrees — MEx (5 26. 

LZPX° / av / av /D>R/RUN/PX°/ av / Hie 72 

ay /D>R/RUN/PZ°/ av / av /D=>R/ 

RUN /+/=/ av / av /R>D/2x° | + |e [29 

“rel 5 30 

Intercept |= CZD —TzD | = [= [31 

(calculated — true zenith distance) | vx [1 [32] 

Post-execution: vy A 33 

/-—/TZD/X/60/=/) arcsin 7 34. 

‘stop 0.35 





(1 in minutes of arc or miles approx.) 
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